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We give the results of density functional calculations for graphene with a widely varying fluorine 
adsorptions. We give a systematic analysis of the adsorption energies, lattice constants, bulk mod- 
ulus, bandgap openings, and magnetic properties. We find that a number of different adsorption 
geometries and a range of physical properties can occur for each adsorbate coverage. The systems 
are found to range from metallic to semiconducting with widely vary band gaps, and a number of 
interesting magnetic phases are found. We expect that many of these structures may occur in real 
materials systems. Further that a listing of the properties found here may help in determining what 
fluorinated graphenes are produced experimentally. 

PACS numbers: 



I. INTRODUCTION 

Since it's experimental realization graphene has been under intense study because of its exceptional physical prop- 
erties including its high electrical conductivity, strength, and elasticity-!"^ However, graphene is relatively inert 
chemically and is a zero band gap semiconductor; thus many of its properties need to be modified for use in ap- 
plications^. Methods to modify the electronic and mechanical properties of graphene include substrate interactions, 
electric and magnetic fields, strain, and chemical functionalization. Surface chemical functionalization is particularly 
a ttra ctive because graphene consists of only surface atoms. A number of adsorbates have been studied, including 
C 8 , and NIT 9 . This work will focus on chemical functionalization with fluorine. 

The fluorination of carbon-based systems has attracted considerable attention. Much work has been directed at 
the fluorination of graphite, particularly in connection with lubrication^. Considerable experimental and theoretical 
work has been done on fluorination of buckyballs and carbon nanotubes^. Fluorine adsorption is a simple chemical 
functionalization that increases the reactivity of sp 2 bonded graphene and opens a wide range of modifications^. In 
addition it enables tuning of band gaps and doping for use in electronics^. Experiments have shown dramatic changes 
in the electronic and structural properties of graphene by increasing the degre e of fluorinatiorf^H^. Theoretical 
work o n fluo rinated graphene has focused mainly on fully covered system s 20 * 21 ! and on systems with low fluorine 
coverage 1 22 * 23 !, with only a few studies that considered several high coverage fluorinated graphene systemiP^HUl 

In this work we consider varying degrees of fluorination (stoichiometrics) on graphene surfaces and a number of 
different arrangements of F atoms (motifs) of the F on surface for each stoichiometry. The overall stoichiometry of 
F on a fluorinated graphene surface is denoted by C m F n where there are n F atoms adsorbed onto m C atoms. In 
Section [IT] we discuss the computational methods used here. In Section |III| the mechanical, electrical, elastic, and 
magnetic properties of fluorinated graphene of varying stoichiometries and motifs are presented. A listing of the 
stoichiometries and motifs studied here and their calculated properties is given in Table [I] We will make frequent 
reference to this table. Concluding remarks are found in Section |TV} 

II. COMPUTIONAL METHODS 

We used the ab initio plane wave code PWscf, which is part of Quantum Espresso, to study supercells of fluorinated 
graphene with periodic boundary conditions®!. The generalized gradient approximation of Perdew-Burke-Ernzerhof 
(PBE) was used to account for exchange-correlation effects. Vanderbilt ultra-soft psuedopotentials^ with wave 
function energy cutoff of 30 R>®1 was employed. We found that an 8-by-8 Monkhorst-Pack mesh sampling of the 
Brillouin zone was sufficiently dense for the total energy to be converged. To speed up the SCF conversion during the 
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FIG. 1: Placement of carbon atoms in a two-by-two unit cell of graphene. The numbers are used to denote the carbon atoms in 
the unit cell that have fluorine atoms attached to them. The dark (light) gray spheres represent the atoms in the A (B) graphene 
sublattice. In the one-sided systems, the attached atoms would be directly over the carbon atoms coming out of the page. 

structural relaxations the wave functions and potentials were extrapolated to second order from the preceding ionic 
steps. 

In this work we studied a large number of periodic fluorinated graphene systems using a 2-by-2 graphene cell as 
the basic unit. This size of supercell is important because as Li et al. have shown, the longest F-F distance that can 
produce a magnetic ground state is about 4.92 A, which corresponds roughly to a 2-by-2 supercelP^l. We will refer 
to the structures studied here by the Item numbers in the first column of Table [I] The second and third columns of 
Table [T] give the stoichiometry and atomic % F of each system. Figure [l] labels the positions of the C atoms in the 
unit cell. The starting configuration of each system had F atoms directly above, or below, the C atoms. The motifs 
of the F atoms in each system are given in the columns labeled "Above" and "Below" in Table [I] where the numbers 
correspond to the C atom having a F adatom attached above or below the graphene planed The calculated values 
of the lattice parameters were obtained by computing the total energy of the F-graphenesystem at a set of lattice 
constants and then fitting the resulting curve with the Birch-Murnaghan equation of stat d 35 * 3 ^ . An inter-layer spacing 
(c lattice parameter) of 24.6 A was used in the direction perpendicular to the graphene planed 

When giving band structures and electronic properties we have taken into account the possibility of ferromagnetic 
and antiferromagnetic spin orderings. When computing the ferromagnetic case, each F atom had an initial spins 
excess of 0.99 hb- In the antiferromagnetic case each of the F atoms adsorbed onto the graphene A sublattice had a 
spin excess of 0.99 fiB and the F atoms adsorbed onto the graphene B had the opposite spirPB The ordering with the 
lowest total energy is the ground state. The adsorption energies, band gaps, intrinsic spin polarizations, and charge 
transfers in Table |T] and in the figures are the values found in the ground state configuration. 

The calculated band structures of all C„F phases studied here can be described by reduced sets of maximally 
localized Wannier functions, MLWF. These are localized, tight-binding-like Wannier orbitals obtained by band in- 
terpolation, which are useful for electronic transport!^ and elasticity models 39 . They describe the occupied bands 
and are a good approximation to the lowest cond uction bands. Wannier90 was used to produce MLWF from the 
Block-like wavefunctions used in Quantum Espresso 40 " 1 !. In performing the Wannicrization, we specified the k-point 
positions of a 4x4 two-dimensional Monkhorst-Pack mesh sampling of the Brillouin zone, which was shifted to include 
the L-point. The band disentanglement was performed over a k-point path from M, to T, to K . With these we can 
reproduce the valence band of C„F and the lowest branches in the conduction band. 

III. RESULTS AND DISCUSSION 

The adsorption energy per F atom versus atomic % F given in Fig. [2] (A) and in Table [I] column 6 may be used in 
understanding the experimental coverage of graphene by F. The adsorption energy per F, Ea was computed using 

E A = {E GR _ nF - (E GR + n F E F ))/n F (1) 

where Ecn_ n F is the ground state total energy of the fluorinated graphene, Eqr is the ground state total energy of 
graphene, E F is the energy of a single F atom, and n F is the number of F atoms in the system. Using this definition 
of the adsorption energyy per F atom the structures with larger adsorption energies have lower ground state total 
energies. We find that the one-sided coverage with the maximum adsorption energy is at 0.2 atomic % F. The dotted 
line in Fig. [2] (A) is used to show the pattern of the maximum adsorption energy for 1-sided motifs. We also see that 
the maximum adsorption energies of the two-sided stoichiometries increase (solid line in Fig. [2] (A)) until full coverage 
(0.5 atomic % F). These two features are consistent with experiment. The one-sided F-graphene motif with the 
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highest per F adsorption energy has F atoms at the para position sites (Table |T] Item 4^21, while the fully fluorinated 
system with the highest adsorption energy is in the chair configuration (Table [i] Item 9l j 24 * 28 l Using modeling of 
the C Ka spectroscopy Bulushcva et al. suggested that C2F is in a configuration similar to Item 26, which we find 
to be one of the most energetically favorable C2F motifs^. Finally, we note that the adsorption energy per F atom 
vs atomic % F is a convenient way of graphically representing other data. We will use it in Figs. [5] and [6] to relate 
conduction and intrinsic polarization properties to the fluorine coverages and adsorption energies. 

We note several points about the calculated absorption energies. For stoichiometrics that have both one-sided and 
two-sided motifs, the motif that is the highest adsorption energy is always two-sided. The one sided (Item 4) and two 
sided (Item 5) C4F coverages with the F in the para positions (adatoms on the C atoms labeled 1 and 6 in Fig. [T]) 
are similar in energy, but when the F atoms are in the ortho (adatoms at 1 and 2) or meta (adatoms at 1 and 7) 
positions the energy difference is largei^. This is likely due to the interactions between the F atoms in the ortho and 
para cases and the relatively smaller interaction in the para configuration. It also appears that the coverages that 
have an even number of F atoms in the 2x2 supercell have higher maximum adsorption energy than their surrounding 
odd coverages. Sublattice imbalance is defined as the absolute difference between the number of F atoms on the A 
sublattice and the on the B sublattice. In the cases where a particular stoichiometry allows for different values of 
sublattice imbalance the motifs that minimize the imbalance have a lower total energy^. 

The optimal lattice constant, as found by the Birch-Murnaghan equation of state, of each motif is given in Table [I] 
column 7 and Figure [2]jB). It shows that the lattice constant tends to increase as the fluorination increases. For 
graphene we found an equilibrium lattice constant of 2.4607 A in good agreement with the experimental value of 
2.4602 A^ForCF, in the chair configuration, we found a lattice constant of 2.61 A which is near the value calculated 
by others^ 7 -* 26 * 27 * 46 !, but is much larger than the experimental value obtained by Nair et al, 2.48 A^. Artyukhov and 
Chernozatonskii claim that this difference can be explained if the surface has a structure like Item 52^. Alternatively, 
this could result from domains of CF of one motif, or different motifs, where the average lattice constant could be 
in the range found experimentally. That domains of different CF motifs could exist is supported by the finding of 
Charlier et al. that a diffusion barrier of about 2.7 eV exists between the chair and boat configurations (Items 51 and 
52), effectively blocking on CF motif from turning into another CF motiP^. Our results also agree remarkably well 
with the C4F lattice constant values reported in Ref. [571 

We also performed variable cell optimizations on a representative subset of motifs (the set of motifs with the 
highest adsorption energy per F atom for each stoichiometry) to determine if the fluorination of graphene caused 
significant changes in the Bravais lattice^. These calculations allow the lattice parameters to be optimized as well 
as the positions of the atoms within the cell. The results gave lattice constants that were approximately equal (to 
within a few hundredths of an Angstrom) to the values given in Table [I] and lattice angles that were consistent with a 
hexagonal lattice. Thus, a simple two-dimensional hexagonal lattice is sufficient for describing fluorinated graphene^. 

When fitting the total energy versus crystal volume curve with the Birch-Murnaghan equation of state one of 
the parameters that can be extracted is the bulk modulus. With fluorinated graphene there is no well defined way 
to account for its thickness. Thus, we cannot accurately compute the bulk modulus. Therefore, Figure [2] (C) and 
Table [I] give the one-dimensional correlation to the bulk modulus as a way to qualitatively discuss the change in the 
elasticity with fluorination. These values were computed by fitting lattice constants versus total energies with the 
Birch-Murnaghan equation of state. Pristine graphene has a elastic modulus of 17.2 eV/A. The values computed for 
fluorinated graphene tend to decrease as fluorination increases. Sahin et al. argue that the rehybridization of the 
carbon atoms into sp 3 structures is the main cause of the change in elasticity 2 ^. 

In Table I] column 9 we list the average C-F bond lengths for each motif (shown in Fig. [2] (D)). Although the 
theoretical lattice constants given in the literat ure for armchair CF (Item 51) are all similar, there is a wide range of 
reported theoretical C-F bond lengths^- 4 * 26 * 28 * 33 !. Our results are near the lower end of the theoretical values reported. 
Our results are in the middle of the range of reported C-F bond lengths for systems like C2F and C4F. Experimental 
values for these t wo sys tems fall over a wide range of values, from 1.36 A to 1.64 A, placing our results within the 
experimental rang 

<fMMM. Covalent (ionic) C-F bonds typically have bond lengths of about 1.32 A (2.01 A). The 
values that we obtain for the fluorinated graphene systems, dc-F = [1-377, 1.554], indicate that the C-F bond is semi- 
ionicP^at low coverage and becomes more covalent as the fluorination increases. This agrees with X-ray photoemission 
spectroscopy whi ch sho ws that the bonds are semi-ionic at coverages less then C2F, but that they become covalent 
at higher coverage 1 51 * 52 !. 

To understand further the nature of the C-F bonds we computed the Lowdin charge transferred from the carbon 
atoms to the F atoms. Table |T] gives the charge transfer per F atom. The charge transfer per F atom versus average 
C-F bond length is given in Fig. [3j There is a general trend of more charge being pulled out of the graphene backbone 
as the adatom bond length increases. We suggest that the increasing charge transfer with increasing bond length 
results from more charge being pulled out of the carbon in order to satisfy the increasing ionic nature of the C-F 
bonds. 

Projecting the DFT results for electronic properties onto tight-binding bands can give physical insight into binding 
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FIG. 2: Calculated values as a function of Fluorine coverage. The gray circles (triangles) represent structures with one-sided 
(two-sided) fluorination. The blue circles (red triangles) are the one-sided (two-sided) motifs that have the highest F adsorption 
energies at a given atomic % F. The values shown are (A) the adsorption energy per F atom, (B) equilibrium lattice constant, 
(C) one- dimensional analog of the bulk modulus (elastic modulus), (D) C-F bond length. The dotted (solid) black line in (A) is 
used to guide the eye in determining the motif of the one-sided (two-sided) structures with the highest adsorption energy per F 
atom. The dashed, horizontal lines in (B), (C), and (D) represent the values computed for pristine graphene. The dotted line 
in (D) is the length of a C-F covalent bond. 



and electrical properties. We do this projection onto a set of maximally localized Wannier functions (MLWF). Figure|4] 
gives the DFT band dispersion plots along with the corresponding MLWF orbital projections of graphene and of several 
illustrative fluorinatcd graphene systems. In the Wannierization of the DFT results the high symmetry energy bands 
are disentangled separately in the following energy windows: the valence band, the localized states in the band gap, 
and the lowest energies of the conduction band. In the valence band, there are branches coming from s-like states 
centered on F atoms that are separated at very low energies (< 10 eV). The other bands are disentangled in terms 
of the following orbitals (Wannier states): one s-like orbital centered on each C-C pair (a bond) including both sp 2 
(non-fluorinated) and sp 3 (fluorinated) carbons, which produce deep energy states; one s-like orbital centered on each 
C-F pair, also with a low energy in the valence band; one p z -like orbital for each pair of neighboring non-fluorinated 
neighboring C atoms (it bond) at the top of the valence band; two p x , p y orbitals on each F atom accounting for 
the paired electrons in the p-shells of the F atoms and giving localized (non-dispersive), low energy states. For the 
conduction bands we interpolate only the branches in a small low-energy window that are relevant for transport 
models. Higher states cannot be disentangled with tight-binding orbitals because they hybridize at high energies with 
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FIG. 3: Charge transferred out of the graphene supercell into a combination of ionic and covalent bonding of the F atoms versus 
C-F bond length. The blue circles (red triangles) represent the one-sided (two-sided) motifs that have the highest adsorption 
energy per F atom for their respective stoichiometrics. 



an unlimited number of quasi-free surface states. We find that the lowest conduction branches are given by p z -like 
orbitals centered on C-F pairs (antibonding a* states) and by a combination of d xz and d^-like orbitals centered 
on unfluorinated C-C pairs (antibonding 7r* states). For fully fluorinated graphene we have only the former set of 
orbitals. 

Figures [1] (A) shows the band dispersion of graphene. The HOMO and LUMO states of graphene are composed 
entirely of the 7r and it* molecular orbitals made from Cp z orbitals. Lower energy valence branches are composed 
entirely of sp 2 molecular orbitals. Due to the symmetry of pristine graphene, there is no hybridization between its ir 
and sp 2 orbitals; therefore its energy bands are not entangled and can be interpolated in a single energy window that 
includes the lowest conduction bands. There are also three high-energy conduction branches created by anti-bonding 
C sp 2 orbitals and a set of quasi-parabolic branches from states outside the graphitic plane. 

Bonds formed between the F p z and the C p z orbitals force the C atoms into sp 3 hybridization, breaking the 
symmetry of graphene. In many cases this opens a gap via the formation of new C-F bonding states deep in the 
valence band and C-F antibonding states in the conduction bancP^ as can be seen in Figs. [4] (B,D,E, and F). 
The resulting bands are entangled due to orbital hybridization, manifested as anticrossings between branches in the 
Brillouin zone. As the degree of fluorination is increased, the number of C-C ir states at the top of the valence band 
is decreased and the number of deep C-F a bonds is increased, until all tt bonds are are removed for CF and the top 
of the valence band is given by a bonds. 

The C4F structure (Item 4) in Fig. [4] (C) is an exception to this trend. The HOMO and LUMO states can be 
completely described by the hybridization of the C pz of the carbon atoms that are not bound to F atoms and the sp 2 
bonding between carbon atoms, though the latter orbitals are only weakly hybridized. The difference in electronic 
structures between the C4F structure discussed here and the other fluorinated graphene structures in this paper 
occurs because the arrangement of F atoms allows for both an aromatic ring and equal numbers of adsorbates on the 
graphitic sublattices. 

For an odd number of non-fluorinated carbons in the supercell the result is a spin-polarized band structure. These 
bands are interpolated by p z -like Wannier states centered on the unfluorinated C atoms. For example, the CsF phase 
has an unpaired p z orbital surrounded by three pairs of n bonds thatgive the valence band in Fig. [I] (B) , and the 
C8F7 has an unpaired p z orbital surrounded by sp 3 carbons as in Fig. Ml (E). 

Table [T] column 11 lists the band gaps of the systems. The letter C is listed instead of a numerical value in cases 
where the system's Fermi level is within a band. Figure [1] (B) shows such a case. Many of the systems have indirect 
band gaps. The systems with indirect gaps will be less photoactive than the systems with direct gaps. 

We find that the majority of the configurations studied here are semiconductors. Figure [5] (A) shows how the 
metallic and semiconducting configurations are distributed as functions of the adsorption energy and the % atomic F. 
The semiconducting systems tend to have larger adsorption energies per F, meaning that it may be hard to produce a 
highly fluorinated system that is metallic. We found that a filled valence band occurred only when the carbon atoms 
without an attached F were found in pairs or rings. Figure [5] (B) shows the distribution of band gaps versus atomic 
% F. We see that it is possible to tune the band gap over a wide range by F adsorption. 

Trends in magnetic properties of fluorinated graphene can be understood in terms of the bipartite lattice of graphene 
and Hund's rule. Each carbon atom in graphene forms three covalent bonds (a molecular orbitals) leaving one electron 
to be part of the tt bonding. If one C atom has a spin up electron in it's p z -orbital then it's neighbors will have spin 
down electrons in their p z -orbitals. The formation of the n conjugation effectively makes the spin polarization of 
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FIG. 4: Band dispersions of: [A] graphene, [B] C 8 F (Table\^Item 1), [C] C A F (Item 4), [D] C4F3 (Item 36), [E] C S F 7 (Item 
44): [F] CF (Item 52). The DFT results are given by light gray, dark gray, and black lines; the black lines represent the spin-up 
results, dark gray the spin-down results, and light gray spin degenerate results. The maximally localized Wannier function 
interpolation is illustrated using blue to represent the C ap i bonds, red the non-bonded C Pz , dark green the symmetric bonds 
between C Pz and F Pz , green the antisymmetric bonds between C Pz and F Pz , and orange the F p bonds. For clarity, [A] and 
[F] are for 1-by-l supercells giving fewer bands than the other figures. 
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FIG. 5: (A) Plot of semiconducting and p-type doped conducting systems. The blue points are for semiconductors and the red 
points are for structures where the Fermi level has band crossings. (B) Distribution of band gaps versus atomic percentage 
fluorination. The blue circles (red triangles) are the one-sided (two-sided) motifs that have the highest adsorption energy per F 
atom. 



the graphene sheet zero. Each neutrally charged F atom will have a spin ±1/2 excess before adsorption on to the 
surface. To satisfy Hund's rule a fluorine's p z electron will have to collapse into a spin state opposite to that of the C 
atom on which it adsorbs. As the F adsorbs onto a C atom it forms a semi-ionic bond, the covalent portion of which 
kills the majority of the spin polarization in the F atom while only slightly changing the spin polarization of the C 
atom. The formation of the C-F bond breaks the 7r-conjugate bonding of the C atom with its neighbors. The nearest 
neigh bors, being the atoms which are most disturbed by the disrupted tt conjugation, share the majority of the spin 
excess^^. The bipartite nature of the surface ensures that the induced spin polarization of the nearest neighbors is 
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shared with the other C atoms that make up the same sublattice^. Further F adsorptions, depending on placement, 
produce strong or weaker induced polarizations in individual atoms. 

Figure[6](A) shows the distribution of ferromagnetic and antiferromagnetic ground states versus F coverage, and (B) 
shows the distribution of motifs that have non-zero net spin polarizations. We find that many of the structures that 
are supposedly in a antiferromagnetic state in fact have a non-zero net spin polarization. Analysis of the intrinsic spin 
of each atom within the different motifs shows a range of magnetic orderings, including nonmagnetic, ferrimagnctic, 
frustrated, and antiferromagnetic. The range of orderings are due to the superpositioning of polarizations discussed 
in the previous paragraph. Items 17, 19, and 20 are an interesting case study. All have four F atoms bonded to C 
atoms 1, 3, 5, and 7, all of the F atoms are on the same side of the sheet in Item 17, while Items 19 and 20 have F 
atoms on both sides. In both Items 17 and 20 the ferromagnetic ordering is the ground state, while Item 17 has a 
strong net intrinsic polarization Item 20 has an almost zero net polarization. Item 19 has three F atoms on one side 
of the sheet and one F atom on the other side. In this case the antiferromagnetic ordering is the ground state with 
an intrinsic polarization of nearly 1 Bohr magnetron per unit cell, which acts like a ferromagnetic ordering. Thus, 
magnetic ordering is highly dependent on the motif formed by the F atoms. 

Figures [6] (C) and (E) show the net spin polarization versus the atomic percent fluorine and the imbalance of F 
adsorbed on each sublattice, respectively. Zhou et al. have calculated the properties of C4H, C2H, and CH, and 
they speculate that, "only half hydrogenation can introduce magnetization" 55 . This is not the case for fluorinated 
graphene. We find that there is at least one motif for every stoichiometry other than fully fluorinated graphene that 
has a nonzero intrinsic spin polarization. Also, contrary to what may be assumed from the simplified explanation of 
intrinsic polarization discussed above, very few of the structures have a net spin polarization that approximately an 
integer multiple of one Bohr magnetron. Although we note that a couple of the C2F motifs are nearly 4 

In order to estimate the thermal stability of the spin polarized states we estimated the polarization energy, Ep, of the 
fluorinated graphene systems (Figurejfj] (D) and Table [i] column 16). The value of Ep is estimated by Ep — E sr — E su , 
where E sr is the energy of the system when computed in a spin restricted calculation, and E su is the energy of the 
system when in the ferromagnetic state (or antiferromagnetic state, whichever is energetically favorable). As is shown 
in the last column of Table [IJ many of the systems have Ep values that are greater than 0.03 eV making it possible 
for them to be stable at room temperature. 

IV. SUMMARY 

We have studied graphene systems for a wide range of fluorinations within density functional theory. Several atomic 
arrangements are possible for most stoichiometrics. A wide range of structural, electronic, elastic, and magnetic 
properties have been calculated for them. Systems vary from metallic to semiconductoring with widely varying band 
gaps. Electronic properties and lattice constants vary with coverage, and a rich range of antiferromegnetic and 
ferromagnetic phases are found. It is reasonable to expect that a variety of these phases may occur in real materials 
systems. 
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